Geometry-Controlled Nonlinear Optical Response of Quantum Graphs 
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We study for the first time the effect of the geometry of quantum wire networks on their nonlinear 
optical properties and show that for some geometries, the first hyperpolarizability is largely enhanced 
and the second hyperpolarizability is always negative or zero. We use a one-electron model with tight 
transverse confinement. In the limit of infinite transverse confinement, the transverse wavefunctions 
drop out of the hyperpolarizabilities, but their residual effects are essential to include in the sum 
rules. The effects of geometry are manifested in the projections of the transition moments of each 
wire segment onto the 2-D lab frame. Numerical optimization of the geometry of a loop leads to 
hyperpolarizabilities that rival the best chromophores. We suggest that a combination of geometry 
and quantum-confinement effects can lead to systems with ultralarge nonlinear response. 
OCIS codes: 190.0190, 160.1245, 020.4900 



I. INTRODUCTION 

The study of nonlinear optical (NLO) properties of ma- 
terials has been the subject of extensive research due to 
the interesting underlying physics as well as the wide 
application range in science and technology. Concerted 
efforts by theorists, synthetic chemists, materials physi- 
cists, and device designers have focused on optimizing 
hyperpolarizabilities, the quantities that govern all light- 
matter interactions, to enable new photonic materials for 
applications in all-optical switching [l[ , 3-D photolithog- 
raphy 0, quantum information, [3[ optical data storage 
[H, photodynamic cancer therapies [f| and other fields. 

In the past decade, the theory of fundamental limits 
found that the first- and second-order hyperpolarizabili- 
ties in the off- resonance regime, @, 0] denoted by f3 and 
7 respectively, depend on the number of electrons con- 
tributing and the energy difference between the ground 
and the first excited states of the system. A compari- 
son of the theory with experiment showed the first hy- 
perpolarizability (which we call the hyperpolarizability 
throughout this text) of the best molecules fell short of 
the fundamental limit (FL) by about a factor of 30. Sim- 
ilar results were found for the second hyperpolarizability, 
7- 

Monte Carlo simulations of the first [§[ and second @ 
hyperpolarizabilities confirm these limits. Based on nu- 
merical optimization of a large variety of potential en- 
ergy functions, the best hy perpol arizabilities are found 
to be about 0.71 of the FL.Jljj, [llj Using a set of exper- 
imental and computational tools including linear spec- 
troscopy, Raman spectroscopy, j3 values measured by 
Hyper-Rayleigh scattering and Stark spectroscopy, Tri- 
pathy et al. studied possible origins of the gap, including 
dilution effects due to vibronic states, unfavorable energy 
spacing, truncation of the sum rules used in the theory 
of fundamental limits and smaller effects such as simpli- 
fications of the Lorentz local field model and systematic 
errors in experiments. Experiments suggest that the na- 
ture of the energy spectrum spacing is the culprit. [T2T] 

The theory of FL predicts that progressively larger en- 
ergy spacing leads to larger /3 values. However, most 



molecules have energy spacing that decreases with in- 
creasing number of states. Recent Monte Carlo simu- 
lations for different classes of energy spectra, including 
hydrogen-like atoms, a harmonic oscillator and a particle 
in a box clearly show that energy spacing is the key fac- 
tor, independent of transition moments, j 1 31] Our studies 
of nanowires are motivated by the fact that a particle in 
the box has favorable energy spacing that is not seen in 
molecules. 

Frohlich first pointed out that at the limit of small 
metal particle size, the continuous electronic conduc- 
tion band becomes discrete. [l4| Kubo used the one- 
electron approximation to show that at low tempera- 
tures, the electronic states of such systems are compa- 
rable to fcT.[lj| Jaklevic et al. reported the observation 
of discrete states in crystalline Pb [l6|, and Brus' model 
of semiconductor crystallites predicted the shape and size 
dependence of the state energies. 17J Ashoori reviews the 
field of artificial atoms. fl8j 

Electrons in quantum wires arc confined in 2-D. Semi- 
conductor nanowires typically have diameters of 1-100 
nm and lengths of several micrometers. Key parame- 
ters, such as chemical composition, diameter, and length 
can be controlled in their synthesis, enabling a wide 
range of devices and applications such as p - n diodes, 
LEDs, transistors and nano-scale lasers. [19j . For exam- 
ple, Huang et al. demonstrated the ultraviolet lasing in 
ZnO nanowires at room temperature, which can be used 
as a nano-scale laser source in optical computing and 
data storage. [2^ Johnson et al. used near-field scanning 
optical microscopy technique to study the NLO proper- 
ties of ZnO nanowires and suggested their use as fre- 
quency converters, pi) Nanowire structures can also be 
used in light generation, propagation, detection, ampli- 
fication and modulation. [22] Yan et al. reported the 
fabrication of programmable architectures using array of 
nanowires. (231] 

Quantum Confined Structures (QCSs) for NLO ap- 
plications were pioneered by Hache et al, who showed 
the enhancement of the second hyperpolarizability of 
gold colloids to originate in the confinement of elec- 
trons and scale inversely proportional to the cube of 
the radius. 24J Sanders and Chang studied the optical 
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properties of Si quantum wires theoretically and showed 
that when the wire's widths is less that 8A, the exci- 
tonic oscillator strength is comparable with the one for 
bulk GaAs semiconductors. [25| In other work, the ef- 
fect of confinement on the luminescence of porous Si was 
emphasized. [2^] Chen et al. showed that for Si quan- 
tum wires, the third order nonlincarity is proportional to 
a -6 , where a is the Bohr radius of the exciton.[27| Other 
studies also confirm the role of confinement in the NLO 
properties of quantum wires. [29j The ability to con- 
trol the confinement in a quantum system enables the 
control of NLO properties of these materials, leading to 
applications in optical devices including lasers, solar cells, 
and nonlinear-optical switches. [30l] 

Past studies have focused only on exploiting the re- 
lation between confinement and nonlinear response of 
QCSs. While the effect geometry on nonlinear response 
of organic materials has been the subject of several 
studies. (3ll - [3^ | no work has addressed this effect on arti- 
ficial systems such as quantum graphs. 

To the best of our knowledge, we investigate for the 
first time the role of geometry on NLO properties of 
quantum graphs, as 2-D networks of quantum wires, 
and show how the hyperpolarizability is affected by their 
shape and orientation. Though past work also consid- 
ered geometry, [35| these were for point charges with 2D 
Coulomb potentials, which are not easily reduced to prac- 
tice as are quantum wires. We show that geometry can 
lead to an enhancement in NLO response. Our present 
focus is on quantum loops as a special class of quantum 
graphs; future work will include more complex graphs. 

We believe we are also the first to investigate the non- 
linear optical properties of quantum graphs and to de- 
velop a consistent, practical method for calculating quan- 
tum observables for graphs. This includes establishing 
methodologies for computing matrix elements of arbi- 
trary quantum operators (including double commutators 
leading to sum rules) between eigenstates of quantum 
graphs. We freely use this lexicography in the text, re- 
ferring the reader to the Appendices for a detailed de- 
scription of the methodology. 

Our paper is organized as follows. First we separate 
geometric and confinement effects by imposing infinite 
transverse confinement to each wire, rendering it one- 
dimensional. Then, we use free-particle one electron 
wavefunctions and energies determined from the loop 
length to calculate the contributions to the nonlinear re- 
sponse of each segment. To test the viability of the re- 
sults, we verify the 2D sum rules for the full loop. Finally, 
we employ numerical simulations to test the configura- 
tion space of geometries to find the largest hypcrpolariz- 
abilitics. We conclude with a discussion of applications 
of this work. 



II. SYSTEM OF A SINGLE QUANTUM WIRE 

A quantum loop is composed of an arbitrary number of 
wire segments. The physics of the quantum graph reflects 
the contributions of single wires. For an extremely thin 
quantum wire, the electron moves freely along the wire 
but it is tightly confined in the transverse direction. This 
tight confinement eliminates size effects .(36j as we briefly 
show below. 

The free-particle wavefunctions along the wire are 



V>(s) = Ae tk ° s +Be 



-ik s S 



(1) 



where s denotes the longitudinal direction; A and B are 
constants; and k s is the longitudinal k- vector which is 
fixed via boundary conditions. 

The uncertainty principle transverse to the wire under 
strong confinement implies infinite energy free particle 
states. [I?) It can be shown (3(| that these states do not 
contribute to the hyperpolarizability by modeling con- 
finement using a highly weighted Delta potential: 



V(t) = - 9 6(t), 



(2) 



where r is the transverse coordinate and g >> 1 is the 
potential strength. The transverse ground state of the 
electron has one bound state wavefunction given by 



(r) = v^e-W, 



where 



, T _ mi 

h 2 ' 



(3) 



(4) 



m being the electron's mass with ground state energy of 
El = ~^. (5) 



2m 



The free state eigenfunctions in the limit of full confine- 
ment, i.e. g — > 00, arc given by 

C en (r) = _yf cosMsin ^ |r A (6) 



and 



€ M {r) = -^~cos(l/7r)sin ^rj , (7) 

where v takes integer values. The even and odd super- 
scripts refer to the wavefunction's parity. Every single 
excited state is doubly degenerate with energy 



E„ 



h 2 k r2 

2m ' 



where kT, is found via 



kla 



VTl. 



(8) 



(9) 
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Here, a is the size of the box transverse to the wire used 
for normalizing the free-particle wavefunctions. When 
evaluating the contributions of the free-particle state, we 
take the infinite space limit a — > oo.[36| 

When a >> 1 Eq. © suggests that the transverse 
excited state energies are positive and become a contin- 
uum, and according to Eq. ([5]), the ground state energy 
approaches — oo. Therefore, the energy gap becomes in- 
finite so these states will not contribute to the nonlinear 
response. [36j 



III. QUANTUM LOOPS 

The typical quantum loop shown in Fig. [T]is composed 
of N wire segments. The wavefunctions of an electron 
moving in this loop are free-particle states with periodic 
boundary conditions and flux conservation along the loop 
[38j . Latin sub or superscripts represent longitudinal and 
Greek characters represent the transverse directions. For 
example E n and E v denote the energy of the n th lon- 
gitudinal and v th transverse states, respectively. Using 
this convention, the Hamiltonian of electron at each wire 
segment is given by 



/V(x) = -g5(xW\ 



.y n+ i) 



(0,0)/ 




(x n , y n ) 



FIG. 1. A 2-D closed quantum loop. (0,0) is the starting 
point of the loop, (x, y) is the electron's coordinate at point v 
with respect to (x n ,y n ). is the angle that this wire segment 
makes with x axis. Each wire segment confines the electron 
in the transverse direction. 



The step- function in Eq. (fT2j) . defined as 



6 [x - x ] 




(17) 



where d x = d/dx and 



d 2 ) - gS (t) = H s + H T (10) 



and 



H T = --d 2 - g S(r) 



(11) 



Using the separation of variables, the solution of the 
Schrodingcr Equation for the Hamiltonian in Eq. (|10[) 

are 

= -J= exp [± (ikfs + #)] 6 [s] 6 [Li - s] ,(12) 



R 

where s is the electron's longitudinal coordinate mea- 
sured from the starting point of the wire segment, 



\J \x - x n f + (y - y n y 



(13) 



R is the total length of the loop, is the longitudinal 
wave number given by, 

, . 2ttZ 



and 



R ' 



= ^ •«((),*- 1), 



(14) 



(15) 



is the accumulated phase from the origin to the begin- 
ning of the wire segment and insures continuity of the 
wavefunction throughout the loop. Si(0,n) is the dis- 
tance along the loop from the origin to the beginning of 



i th wire segment. 

i-l 



s(0,i- 1) = J2 L i = H V ( x o - ^-i) 2 + (Vi ~ %-i) 2 - 



3=1 



3=1 



(16) 



demands that the wavefunction on i th segment vanishes 
everywhere except along the i th segment. 

The longitudinal energy of the electron is found using 
periodic boundary conditions, ip(s + R) = tp(s), yielding 



2n 2 h 2 l 2 
mR 2 



Z = ±l,±2,--- 



(18) 



Given the symmetry between clockwise and counterclock- 
wise motion, the longitudinal wavefunction is doubly de- 
generate. 

The transverse wavefunctions of the electron at each 
wire segment are given by Eqs. (j3J), © and ([7]). The 
total wavefunction of the electron is the union of all seg- 
ments including the transverse and longitudinal parts. 
Using the notation outlined in Appendix \K\ the total 
wavefunction, </> m ^(s,T), can be written as a union over 
the graph of the wavefunctions 0^ m (s,t) on each edge: 



</w(s,t) = Vm(s)^(r) 



(19) 



The total wavefunctions are orthogonal because they 
span a Hilbert space defined by their Hamiltonian. How- 
ever, the value of the total wavefunction on a particular 
segment has no quantum mechanical significance by it- 
self: All such contributions arc required (with matched 
boundary conditions) to obtain the energy spectrum from 
a secular equation. Consequently, the value of an eigen- 
fucntion on a particular segment depends on its value on 
all other segments. Two different modes on a particu- 
lar segment don't have any special relationship between 
them. 

Each energy eigenvalue is a sum of longitudinal and 
transverse energy components, given by Eqs. (|18|) and ([5]) 
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for the transverse ground state and |5J for the transverse 
excited states, respectively. 

The longitudinal excited state wavcfunction are dou- 
bly degenerate because of the clockwise/counterclockwise 
symmetry of the Hamiltonian, and the transverse excited 
state wavefunctions are doubly degenerate and are eigen- 
states of the parity operator, with odd and even eigen- 
values, in Eqs. © and (J7J, respectively. 



IV. TRANSITION MOMENTS 



The hyperpolarizability is a function of energies and 
position matrix elements, x^'s, which we loosely call 
transition moments. The transition moments of an elec- 
tron in a quantum loop are related to the sum over con- 
tributions from each segment by (see Eq. (|T9| ) 



X P K,q\ = (pK\x\qX) 



(20) 



dxxcff (x)(j) q x{x) 



E 

k 



dxxi 



b k p l(x)4> k qX {x). 



For the i wire segment, the x coordinate of the elec- 
tron with respect to the starting point of the graph, (0, 0) 
in Fig. [TJ is given by 



s cos 9i ± r sin ( 



(21) 



where 6i is the angle that segment i makes with the x 
axes, and x\ is the x coordinate of starting point of the 
segment i. '±' represents the two possible sides of the 
wire from which r can be defined. The total transition 
moment is given by substituting for Xk in Eq. (|20|) and is 
the sum of a longitudinal contribution and a transverse 
contribution. 



Eqs. (|22j) and (|23|) suffice to calculate x^-'s for any wire 
orientations including when x\ = X2 or yi = y 2 - 



2. Transverse Transition Moments 

We follow the same procedure to find the contributions 
from the transverse component. Due to parity selection 
rules, only odd (even) to even (odd) transitions are per- 
mitted. At each wire segment, 



4) 



pn,q\ 



sin6>i / dsipUs)ip l Js) / d-q-q-q K (T)rj\(T) 



= sin0; x P pq x J l KX 



where 



Lx 
R 



ik 



(0,i-l) / 



p = q 
1) P^Q 



(24) 



(25) 



and 



r -8V2(-i)^^/^3 /2 >o a(k) = 

(a 2 fe; 2 +47r 2 K 2 ) V ; ' \ J 



2(l + (-i; 



a | k. A | 



7T 2 (K 2 -A 2 ) 2 



|«| = |A|^0 



(26) 



The graph's transition moments, Eq. (|2"Tj) , are ob- 
tained by summing over longitudinal and transverse tran- 
sition moments from all wire segments, given by Eqs. 
m or ([231), and ([24]). 



V. SUM RULES FOR QUANTUM LOOPS 



1. Longitudinal Transition Moments 
For the two distinct cases given by p =/= q and p = q, 

{4) pK , qX , Eq. m Yields, 

a: p^q: 



( x i\ = £^ e « w -(o,<-i) x 

\ s /pK,q\ Jj> 



cos (9; in \ -t, j cos 9i x\ 

Kqp J Kqp K,qp 



U2 



Note that i 
the index i. 



(22) 



T should not be confused with 



b: p = q: 



(4) 



\x\ + x\J Li 
2 ~R 



(23) 



The Thomas- Reiche-Kuhn sum rules are a consequence 
of the commutators between the position operator and 
the Hamiltonian of the quantum system, [x, [H,x]], and 
relate the transition moments and energies. In 1-D, the 
sum rules are given by 



Em. + E r , 



2m 



(27) 



In general, the summation spans the complete set of 
cigenstates of the system, including both discrete and 
continuum [3^| and degenerate and non-degenerate states 

^ i-l 

Sum rules have been widely applied. [41[ In nonlinear 

optics, they are used to calculate the fundamental lim- 
its of first and second hyperpolarizabilities,[(| Q to find 
a dipole-free sum over states expression for calculating 
the nonlinear hvperpolarizabilities (42l l43l | . and to impose 
constraints on the related parameters in Monte Carlo 
simulations of the first Q and second [j| hyperpolarizabil- 
ities. In this section we obtain an analytical relation for 
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TABLE I. Numerical results for diagonal components, i.e. p = 
q and k = A in Eq. (|28fl for triangles, 100 longitudinal and 
200 transverse states. 



Triangle Coordinates 
(x ,ya,xi,yi,X2,y2) 


Longitudinal Transverse 
(xh 2 /m) (xh 2 /m). 


(23,-16,3,-19,24,-15) 

(18,-7,21,1,-25,-20) 

(-32,-21,-49,-42,-69,-88) 

(3,-93,-13,36,35,20) 

(-7,-16,-27,64,42,-23) 


0.474 0.022 
0.401 0.095 
0.118 0.378 
0.093 0.403 
0.196 0.299 



sum rules in the longitudinal and transverse directions of 
a quantum graph and use these to validate the calculated 
eigenfunctions. 

Sum rules in 3-D were verified for the quantum rigid ro- 
tators by Hadjimichael et ai, who showed that the classi- 
cal picture of a rigid rotator, where the radius of rotation 
is fixed, is consistent with the sum rules when taking into 
account the radial component of the wavefunction. [37| 
They used the uncertainty principle in the radial com- 
ponent under a highly-confining delta function potential, 
V(r) = —gS(r — R ) with g > and R being the po- 
tential strength and radius of the rotator, respectively, 
to show that the continuum states corresponding to high 
radial momentum ensure that the sum rules are obeyed. 

In reference [36| , it was shown that Eq. (f27|) for a single 
quantum wire may be written as a sum rule over s plus 
a sum rule over r. Since the sum rule is derived from a 
double commutator of x with H, it is an operator which 
may be evaluated across the graph using Eq. (|A10I) in 
Appendix [X] as a sum over edges of the sum rules for each 
wire segment: 



E 



El 



2 J ^pn^nq 



- 2 



E r _ K+ E X \ ri r,i T,i 
2 / VI KV v\ 



(28) 



h X 
2m 



where P pq is the edge overlap integral defined in Eqn (f25|) . 
In deriving Eq. (|28]l. we used the closure relation Eq. 
(Hl|to write EnjVl^Ipr 

The first and second terms on the left-hand side of Eq. 
(|2"8")) can be interpreted as the longitudinal and transverse 
contributions, respectively. 

We can use Eq. (f2"5| to test the validity of the re- 
sults in section IIII1 We pick the simplest possible loops, 
i.e. triangles, made of three quantum wires connected 
at vertices. Many configurations are tested by randomly 
sampling the vertices. The energy and transition mo- 
ments for each are calculated as previously described, 
and evaluated by the expression on the right-hand side 
of Eq. (|28p . The numerical result for six different trian- 
gles are presented in table U for diagonal sum rules only. 



Verifying non-diagonal sum rules follows the same proce- 
dure. The sum of second and third columns in Table HI 
when multiplied by h 2 /m yields the left-hand side of Eq. 
([2"8")l . thus confirming the validity of our protocol. 

The projections of the longitudinal transition moments 
on the lab coordinate axes depend on the geometry 
through the orientations of the wire segments. Conse- 
quently, the nonlinear optical response of quantum loop 
depends on the shape and its orientation. In contrast, the 
projections of the transition moments due to the trans- 
verse and longitudinal coordinates together leave the sum 
rules invariant to the loop's shape or orientation. 

The transverse wavefunctions do not contribute to (3 
and 7 because of the infinite energy denominators. How- 
ever, they make finite contributions to sum rules because 
the product of infinitesimal transition moments and in- 
finite energy EijXij is finite. Hence, when calculating 
the hyperpolarizability we only use the longitudinal parts 
given by Eqs. (j2"2"|) and 



VI. 



GEOMETRY-CONTROLLED NONLINEAR 
RESPONSE 

A. First Hyperpolarizability 



In this section, we investigate the effect of geometry 
on the hyperpolarizability of a quantum graph. The off- 
resonance diagonal component of the hyperpolarizability, 
P xxx , is given by 0], 



(29) 



where x\j is the transition moment of the quantum graph 
along the coordinate axis x. x nm = x nm — xooS nm and 
the prime in the summation indicates that the ground 
state is excluded from the summation. The fundamental 
limit formalism defines j3 to be the largest diagonal tensor 
element. For a given quantum system, /? is determined by 
finding the coordinate frame in which j3 xxx is the largest. 
In the present work, we do so by fixing the coordinate 
axes and determining /3 XXX as a function of angle, from 
which we determine the peak value /?. One can take the 
same approach to calculate the effects of geometry on 
other components such as /3 xxy and f3 xyy . 

For a quantum system with N electrons and energy 
difference between first excited state and ground state, 
Eio, the theory of FL predicts that the largest attain- 
able diagonal element of the hyperpolarizability, (3 max , is 
given by Q 



,1/4 



eh 
ZJ/2 



N 3 / 2 

p 7/2 



(30) 



To optimize the nonlinear response regardless of the 
molecular size, we use the intrinsic hyperpolarizability, 
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P™ x = -0.049 







-3 

X 



FIG. 2. The three triangles with largest /3 int (~ -0.049). 



Pint: as a scale-invariant figure of merit, [4 



3\ X ^'^OnCranCmO 



Pmax V 4 7 e 6 n e m 
x ' n.m 



(31) 



where ^ and e, are normalized transition moments and 
energies, denned by 



^ ~,rnax 5 
A 01 



where 



' in 



\2mE w ) ' 



(32) 



(33) 



and where xJJ| a:c represents the largest possible transition 
moment, xqi. According to Eq. (|32|) . eo = and ei = 
1. /3i„t is scale- invariant and can be used to compare 
molecules of different shapes and sizes. 

We sample triangle configurations by randomly choos- 
ing the vertex coordinates, then calculating the total en- 
ergy and transition moments using Eqs. (|T5|) . (|2"5|) and 
(|22p. Subsequently, /3, nt is calculated using Eqs. (j3"Tj) and 
(|32[). For each configuration, the triangle is rotated in x-y 
plane in increments of 15° and (3 int is calculated for each 
orientation. The largest Pint value and the correspond- 
ing vertex coordinates are recorded. After thousands 
of trials, the best intrinsic hypcrpolarizability found is 
\P int \ = 0.049. 

Figure [5] shows the three triangles with largest Pint — 
—0.049. Though the triangles differ in area, the scale 
invariance of Pi„t assures that they can be compared. 
Rotating each triangle by 180° identifies the orientation 
with peak (3i nt of positive sign. 

Figure [3] plots (3 XX x versus the ratio of the smallest to 
the largest internal angles, m i n /8max - a measure of the 
shape of a triangle. The blue points label the isosceles 
triangles. 

Since simulations show isosceles triangles with the 
largest hyperpolarizabilities, we more systematically 
study them by fixing two of the vertices at coordinates 
(0,10) and (0,-10) while the third vertex is scanned 
along y = for x 3 > 0. Fig. [4] shows a plot of pint as a 
function of x% yielding a peak of Pint ~ 0.049 for the tri- 
angle in the inset. This peak value is the same for all tri- 
angles with two angles ~ 68° and ~ 44°, independent of 



ca 



0.050 



0.025 



0.000 



-0.025 - 



-0.050 - 




e . .e 

mm/ max 



FIG. 3. Distribution of Pint over its full range of triangular 
configurations as quantified by the ratio of the smallest an- 
gle, Omin, and largest angle, 9 max . Points in blue represent 
isosceles triangles. 



ca 



0.05 


i i i i 




0.04 




44^» " 


0.03 






0.02 






0.01 






0.00 


i i i i 





20 



40 



60 



80 



100 



FIG. 4. Pi n t for isosceles triangles with two points fixed at 
(0,10) and (0,-10), and the third point varying along the x 
axis. The triangle in the picture represents the largest value 
of p int = 0.049 



their size. Previous studies that varied the positions and 
charges of 4 point nuclei found that the optimum geome- 
try with 72° and 35° yielded Pint « 0.65. [35| However, it 
is not possible to control point nuclei with such precision, 
let alone synthesis molecules with arbitrarily charged nu- 
clei and bond angles. This is why real molecules fall far 
short of the limits. 

Studies which optimize the shape of the potential en- 
ergy function yield the largest possible nonlinear optical 
response of Pi n t — 0.71 for a large set of potentials, which 
universally have the property that X = xqi/x™^ where 
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x rnax j g gj ven by JTJq 



©•(IS El El Since the opti- 
mized triangle has X = 0.49, isosceles triangles are not 
globally optimized, but represent a local maximum un- 
der constraints of their geometry. In future studies, we 
will consider other graphs, including bent wires, three- 
pronged star graphs and fractals and will show that a 
much larger value of j3i n t possible. 

In the three-level ansatz, p int can be expressed as.pTj 



o.ooo 



where 



Pint = f(E)G{X), 



f{E) = (1-Ef/ 2 [E 2 + ^E+l 



G(X) = VzX x r-{l-X% 



(34) 



(35) 



(36) 



and E = E-io/E^o- The energy spectrum for a trian- 
gle is similar to that of a particle in a box, which has 
large intrinsic hyperpolarizability. [l3| For a triangle with 
E = 3/8 and X = 0.49, in principle, the hyperpolariz- 
ability could be as high as Pint = 0.65 for a three-level 
system. Comparing this value with p in t ~ 0.05 suggests 
that a larger number of states contribute to the nonlinear 
optical response. Also, FL theory predicts that degener- 
acy lowers the response. /3™' for the optimized triangle 
shows that 5 states contribute significantly. These two 
factors suppress the full potential of a triangle. 

The search for best NLO quantum loops can be gener- 
alized to 3+ segments. As for triangles, the coordinates 
of four vertices are chosen randomly without restrictions. 
In configurations where wire segments cross, we assume 
that there is no transfer of probability current between 
the two - for example one segment can be viewed as be- 
ing above the other one by an infinitesimal distance with 
tunneling effects neglected. The best quadrilateral gives 
Pint = 0.073, which is 50% larger than for the best trian- 
gles. For a larger number of segments, the improvement 
in Pi n t is marginal. Future work will address more com- 
plex graphs including many-electron band-filling models. 



B. Second Hyperpolarizability 

The sampling approach for calculating p for loops may 
be applied to the calculation of the second hyperpolariz- 
ability, 7, using the same matrix elements that we used 
for p. In these studies, we seek to optimize the intrinsic 
second hyperpolarizability given by 



lint 



4 I £—1 p_p__p, £-~< 



l n.m.l n.m 



n t-m 



(37) 




-0.150 



where the normalized transition moments, and ener- 
gies, ei, were introduced in Eq. I32I 



e . ,e 

mm/ max 



FIG. 5. Distribution of 7i„t over its full range of triangular 
configurations as quantified by the ratio of the smallest to 
largest angle. 



The result for a triangle loop yields a maximum value 
Ixxxx = for high aspect ratio (flat or squashed) trian- 
gles lying mainly along the y-axis and the largest negative 
value of -0.138 for similarly-shaped triangles lying nearly 
parallel to the x-axis. Fig [5] displays the variation of 
Ixxxx over the entire range of aspect ratios for all 10,000 
triangles sampled in the Monte Carlo run. 

The two branches correspond to the minimum and 
maximum values of the intrinsic hyperpolarizability for 
a given aspect ratio as quantified by m i n /0 max . For an 
equilateral triangle, m i n /6 rna x = 1, the two branches 
merge to yield the same minimum and maximum val- 
ues simultaneous. At the other extremum of a large as- 
pect ratio triangle, the minimum and maximum values 
of ixxxx are at the extremes, suggesting that large as- 
pect ratio leads to the largest magnitude of the second 
hyperpolarizability. 

Similar results hold for loops with four edges. The sec- 
ond hyperpolarizability is always negative, and its min- 
imum is at about -0.138. This remarkable result illus- 
trates that graphs with identical topologies but different 
geometries have similar nonlinear optical responses. This 
is the topic of a future paper. 

Comparison between Figures [3] and [5] reveals that there 
may be a relationship between the two. In both branches, 
the magnitude of Pi„t increases as the triangle becomes 
open. In comparison, the largest magnitude of jint de- 
creases as the triangle becomes more open. Thus, there 
may be a complementarity between the two, as is often 
found in other properties. 



VII. CONCLUDING REMARKS 

Low dimensional quantum systems are scientifically in- 
teresting because of their novel quantum mechanical ef- 
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fects, and should have a broad range of potential appli- 
cations in future technologies. The fact that the optical 
properties of these systems can be manipulated by vary- 
ing the degree of confinement provides a useful tool for 
the fabrication of novel nano-scale structures for new de- 
vice concepts and to provide novel approaches for making 
more efficient nonlinear materials. The present work has 
established a theoretical foundation for using quantum 
graphs as prototypes for investigating how geometry im- 
pacts the NLO response. 

We deliberately prepared a single quantum wire in a 
way that made confinement effects irrelevant to the non- 
linear response. Then we solved the Schrbdinger Equa- 
tion for an electron in a wire as part of a larger net- 
work. The 2-D sum rules were obtained and were used 
to test the validity of the results. It was shown that 
the transverse wavefunction of the electron is essential 
in validating sum rules, but plays no role in the NLO 
response. Thus the hypcrpolarizability of the quantum 
graphs with different shapes and identical topology are 
solely controlled by geometrical effects. 

Monte Carlo simulations identified the best triangle 
loops with the largest /3i nt values as a unique class of 
isosceles triangles with orientations and internal angles 
given in Fig. 01 4-scgment loops show even larger 
values. Preliminary results for rectangular loops show 
an intrinsic value near 0.07. While this value is far 
from the fundamental limit, this geometry-controlled 
value is better than the hyperpolarizability of all other 
known nanosystems. This study, which focuses on ge- 
ometrical determinants of [3i„ t and 7,„t shows promis- 
ing results toward replacing nonlinear optical materi- 
als with artificially structured systems. Recent work on 
more complex graphs have yielded much higher intrinsic 
hyperpolarizabilities. [HI 

We also showed that the triangle graphs have negative 
jint- This appears to be a consequence of the closed- 
loop topology of the graph and not its geometry. Parallel 
studies underway in our group have shown that opening 
a vertex of the triangle without changing its geometry 
allows for positive values of j XX xx- 

The next potential direction in exploring the nonlinear 
optical response of quantum graphs would be generaliz- 
ing the one-electron model to N-electrons, and consider- 
ing the resonant regime, where the theory can be com- 
pared with experimental results and be more relevant 
for practical applications. Meanwhile, it is also interest- 
ing to study more complex graphs, from the more exotic 
three-pronged star graph [49| to fractal systems, where 
the extra degrees of freedom lead to greater design flexi- 
bility. 

Additionally, the mathematical formalism that we have 
developed will be of value in evaluating more complex 
graphs by relating the properties of individual nanowircs 
to the properties of complex structures in which they are 
fundamental building blocks. 
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Appendix A: Quantum graph lexicography 

A quantum graph is a metric graph supporting particle 
dynamics. It is characterized by a Hamiltonian H with a 
complete set of eigenstates \N) and an energy spectrum 
£jv satisfying H\N) = E^\N). If the graph is planar 
and the particle is confined to move on the E edges, the 
graph is a quasi 1-D system. We fix the graph in the 
xy-plane. 

The energy eigenstates \N) span a Hilbcrt space and 
are constructed from a union of non-intersecting sub- 
spaces, each spanned by vectors \N k ), where k labels an 
edge, which also satisfy H\N k ) = E N \N k ). The eigen- 
states \N) may be written as a union of edge states with 
the interpretation that the dynamics of the particle on 
edge k is described by the set of edge states \N k ) whose 
projections into real space (x\N k ) = (^(x) are wavefunc- 
tions exactly equal to the projection of the eigenstates 
4>n(x) = (x\N) on edge k. The eigenstates and oper- 
ators acting on them comprise the Hilbert space of the 
graph. 

Graph dynamics, including nonlinear optical phenom- 
ena, are calculated from the eigenstates. Quantum me- 
chanical relations, such as matrix elements or sum rules, 
are expressed in terms of eigenstates, not edge states. 
Edge states are not complete on the graph: A set of 
edge states on the k th edge represents dynamics on that 
specific edge, not on other edges. This implies that in- 
ner products of eigenstates will necessarily be sums of 
inner products of edge states from identical edges. In 
real space, this simply means that the integral over the 
graph of a bilinear form of eigenstates with any opera- 
tor is parsed into a sum of integrals over edges of bilin- 
ear forms of edge wavefunctions. Bilinear forms such as 
(A l |M J ) are meaningless unless i = j. 

The edge wavefunctions form a complete set on an edge 
only when every eigenfunction on an edge obeys the same 
boundary conditions at the ends of the edge. Eigcnfunc- 
tions that do not obey such boundary conditions are not 
complete on an edge, but, the wavefunction composed of 
the union of all edges produces a complete set of eigen- 
functions of the Hamiltonian on the full graph. The liter- 
ature often refers to edge spaces as Hilbcrt spaces, though 
they are actually subspaces that only in union produce a 
Hilbert space representing a dynamical quantum system. 

In position space, solving a graph means finding 
the edge wavefunctions and energy spectrum, matching 
boundary conditions at the vertices, and constructing the 
cigenfunctions as a union of edge wavefunctions using an 
appropriate notation. This approach ensures that the ab- 
stract definition of the quantum graph is realized prop- 
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erly in position space. All physical quantities may be 
calculated using the edge wavefunctions and the energy 
spectrum, provided that the definition of the union op- 
eration is properly implemented. 

To this end, it is useful to define the union operation 
in abstract space and use it to derive expressions for ma- 
trix elements of operators, inner and outer products of 
eigenstates, and sum rules. All of these may be reduced 
to sums over the graph of integrals of bilinear products 
of edge wavefunctions and real-space operators. But the 
abstract form enables almost instant proof of theorems 
such as sum rules for graphs for which proof has been 
shown for single-edged graphs (quantum wires). 

We emphasize here that the following abstract defini- 
tions are exactly that - abstractions defined to facilitate 
the manipulation of edge states comprising true eigen- 
states when physical properties are to be calculated. The 
abstraction is merely book-keeping; the definitions could 
easily well be written using alternate notation. 

We define an energy eigenstate vector as 



|A0=uf =1 |A fe ) 



and the dual (bra) state as 



(M| = (M^|u]f l5 



(Al) 



(A2) 



where the union is over all the edges, E. The operation 
U*f = i i s placed to the left of the ket while its conjugate 
is placed to the right. We do this for convenience and as 
a reminder of what gets 'unioned' to what. The union 
operation is not an operator in any sense ordinarily used 
in Hilbcrt space; it is an instruction that can be placed 
anywhere that is convenient, and can be interpreted ac- 
cordingly. 

The union operator always accompanies an edge bra 
or ket in any expression involving more than one edge 
of the graph. Note that the Hamiltonian operating on 
an energy eigenvector results in H \N) = U^ =1 H \N k ) = 
En Uf =1 \N k ) = E N \N) since all edges have the same 
energy spectrum as the graph. 

The inner product of two energy eigenstates in this 
notation is then 



(N\M) = (M'| utf, Uf =1 \N l ) = 6 NM 



(A3) 



since the energy eigenstates are orthonormal. The inner 
product is defined as: 

(N\M)= J dsiP*(s)^(s) = Y, J dx k ^ {x k )4> k M {x k ) 



E (M k \N k 



k=l 



(A4) 



Equations I A3I and I A4I define the allowed operations from 
which all others follow. 

Similarly, we may compute the outer product of two 
eigenstates as 



\N)(M\ 



uf^lM^lutfi 



(A5) 



When Eq. (|A5|) is inserted between two operators or 
eigenstates, the union operations act to select out the 
proper edges on either side of the outer product to ensure 
that all physical quantities are calculated on the same 
edge, then summed over edges. 

The definitions in Eqs. (|A1|) and (|A2I) may now be 
used to express the matrix elements of an operator as 
sums over edge states in a way similar to the derivation 
of Equation 



<JV|0|Af> = (N j \ utf x Ouf =1 | AT) 

E 

= ^2(N*\0\]Vr) (A6) 



The completeness relationship in Eq. (|A5|) may be 
used with the definitions above to derive the matrix ele- 
ment of a composite operator AB by inserting it between 
AB. It is equivalent to the standard approach of inserting 
a complete set of eigenstates between AB and using the 
definitions of the bra and ket expansions in edge states 
to find: 



(N\AB\M) 
J2(N\A\P)(P\B\M) 



(A7) 



E (^'1 | J*) (P k \ U& BujLi \M l ) 

p 

E \J2( Nj \ A \ pj )j (x><ii?iA/^ 



where the last step follows from Eq. (|A3|) and the fact 
that the unions do not act upon the operators themselves. 
This approach to expressing matrix elements of compos- 
ite operators in sums over edge state matrix elements 
may be scaled repeatedly to any number of operators. 

To compare the use of eigenstates and edge states, con- 
sider the derivation of the the sum rules. The sum rules 
follow from calculating the matrix elements of the com- 
mutator [[H. x], x].(50t 



(N\ [[H,x],x] \M) = (N\Hxx + xxH ~ 2xHx\M) 



2m 



>NM ■ 



(A8) 



Using closure, the first term may be written as 
(N\Hxx\M) =J2E N (N\x\P)(P\x\M). 



(A9) 



Expressing the other terms in this same way and sum- 
ming them leads to the usual TRK sum rules. [H3|. 

It is now a simple exercise to express the TRK sum 
rules in terms of edge states by using Eq. (|A6[) with 
O = [[H,x],x]. By using the closure relationship Eq. 
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5|) between pairs of operators, the result is that 



Appendix B: Calculation of transition moments 



E ( - U E » + Em)](N 1 \x\P*)(P3\x\M1) 



1:3 = 1 \ P 

fe2 E 



2m 



(AlO) 



It has been shown that the TRK sum rule decomposes 
into a long itudinal and transverse part on a single wire 
(edge) |36| by using the commutator of H with both lon- 
gitudinal and transverse coordinates for the wire. The 
general term on the left hand side of Eq. (|A10[) may be 
written as the product of matrix elements between edge 
states of the double commutator, since edge states are 
themselves eigenvectors of H. Thus, the decomposition 
of a single wire sum rule into longitudinal and transverse 
parts holds exactly for the sum over edges of all such 
wires. 

This formalism can be generalized to more complex 
graphs and curved edges by either converting sums to 
integrals in the infinitesimal edge limit, or by defining a 
curvature metric that maps an edge to a curve. 

A final remark about edge functions: They must not 
be used on their own to compute global quantities for 
a graph, but must always appear in unions over edges. 
However, it is likely that graphs with similar motifs em- 
bedded in their topology will have edges with similar 
wavefunctions. As an example, a graph with two loops 
separated by a single wire (a so-called barbell graph) 
has edge wavefunctions on the loops that are identical 
in form and differ only in details such as edge lengths. 
When writing a set of solutions from scratch for a new 
graph, one can write down the form of the edge functions 
almost immediately and even select their amplitudes so 
they match easily to the connecting wire. The specific 
relations among the amplitudes for the loops and con- 
necting wire are then found from the conservation of flux 
equations, which also lead to the energy spectrum for the 
graph. It is unknown at this time whether or not there 
are sum rules associated with individual edges. 



As represented in the above formalism, the absolute 
position and orientation of an edge relative to the lab 
frame is not explicitly labeled in an edge kct in E-spacc. 
In the an edge's frame, we define s along the edge and 
f perpendicular to s (i.e. Eq. (TIT?]) with — > s and 



y 



f 



t). We limit ourselves to two dimensions, though 
edges in higher dimensions are treatable in the same way. 
The wavefunction along an edge k is (f>k{s) = (s\ (j) k ). 
The same holds for r. We may choose either end as the 
starting point. 

When calculating the expectation value of an observ- 
able in the lab frame, we can transform the wave function 
to the lab frame, or we can transform the observable into 
the edge frame. We choose the latter. For edge k with 
the reference end at x\ and y\ , the position operator can 
be expressed as, 

x k = tx k + l T s cos 6 — t s T sin 9 

y k = tx\ + l r ,ssm6» + 1 s tcos6>, (Bl) 

where t s and t T are the identity operators in the sub- 
space s and t. 

The contribution from edge k to the position matrix 
element is given by our conventions as defined above, 



( P k\ x k \qX) = 5 kA / dstp k *(s) (x k + scos9) ^ k (s) 




ih^(r), (B2) 



where p and q are quantum numbers for states along the 
edge; k and A are quantum numbers for the transverse 
states, and L k is the length of edge k. The edge overlap 
integral I k q was defined in Eq. (|24|) . The second term 
was shown to vanish for tight transverse confinement. (36j 
This illustrates the derivations of Eqs. (|2"2"]l and (|23p . 
with the definition {xl) pKqX = {pn\x k \q\) and a mod- 
icum of algebra. 
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